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Preface 

In  this  work  I  have  Bade  an  effort  to  model  the  complex 
power-burning,  chambered-gun  internal  ballistics  process  with 
an  artificial  sequence  of  fundaaental  processes.  Although  I 
was  unable  to  obtain  acceptable  results  from  the  aodel  in  the 
time  allotted  for  this  work,  I  feel  that  the  partial  results 
attained  to  date  indicate  that  the  aodel  shows  good  promise. 

At  the  very  least  I  have  learned  a  great  deal  about  the  gun 
business,  the  application  of  engineering  principles  to 
physical  problems,  the  value  and  results  of  simplifying 
assumptions,  and  the  frustrations  involved  in  creating  and 
perfecting  a  lengthy  and  involved  computer  program. 

I  now  take  this  opportunity  to  express  ay  gratitude  to 
my  thesis  advisor.  Dr.  James  Hitchcock,  both  for  posing  this 
most  challenging  problem  as  well  as  for  his  knovledgable 
advice  on  analyzing  the  gun  problem.  I  am  also  grateful  for 
the  timely  suggestions  m^de  by  the  other  gentlemen  on  my 
thesis  committee.  Dr.  Andrew  Shine  and  Capt  (Dr.)  Stephen 
Koob.  And  I  thank  my  lovely  wife  Judy,  whose  patient  under¬ 
standing  during  this  difficult  time  has  been  truly  remarkable. 

James  F.  Setchell 
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Abstract 

An  attempt  is  made  to  model  the  internal  ballistics 
process  of  a  powder-burning  gun  by  replacing  the  actual 
internal  ballistics  process  with  an  incremental  sequence  of 
phases.  These  phases  are  a  constant-volume  energy  transfer 
phase,  a  shell  motion  and  finite-amplitude  wave  propagation 
phase,  a  propellant  motion  phase,  and  a  gas  expansion  and 
mass  transfer  phase,  The  model  permits  consideration  of  a 
chambered,  powder-burning  gun  problem  with  unspecified 
pressure,  density,  velocity,  temperature,  and  propellant 
distributions.  The  method  of  solution  shows  promise,  but 
useful  results  have  not  been  attained  to  date. 
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AN  ATTEMPT  TO  MODEL  THE  GUN 
INTERNAL  BALLISTICS  PROBLEM 

I .  Introduction 

Background 

Gun  internal  ballistics  is  the  study  of  the  conversion 
of  latent  chemical  energy  of  a  propellant  to  kinetic  energy 
of  a  projectile.  It  is  only  concerned  with  the  period  of 
time  that  begins  with  propellant  ignition  and  ends  with  the 
projectile  leaving  the  barrel.  The  primary  purpose  of  a 
gun  internal  ballistics  study  is  to  predict  the  gas  property 
and  shell  motion  history  of  a  gun. 

The  formal  study  of  gun  internal  ballistics  began  during 
the  eighteenth  century  with  the  work  of  Benjamin  Robins  in 
1742  and  Count  Joseph-Louis  Lagrange  in  1793.  C.  K.  Thornhill 
includes  a  summary  of  early  gun  internal  ballistics  work  in, 

"A  New  Special  Solution  to  the  Complete  Problem  of  the 
Internal  Ballistics  of  Guns",  and  suggests  that,  since  the 
time  of  Robins  and  Lagrange,  analyses  of  the  gun  problem 
have  generally  followed  one  of  three  methods.  The  first 
method  involves  a  solution  to  the  complete  fluid  dynamic 
equations  of  flow  using  the  theory  of  finite-amplitude  waves 
in  gases.  The  second  method  involves  a  reduction  of  the 
problem  to  the  solution  of  ordinary  differential  equations. 
Solutions  employing  this  second  method  are  known  as  "mixed 
solutions",  and  do  not  involve  the  complete  fluid  dynamic 
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equations  of  flow.  Solutions  Involving  the  third  method  are 
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known  as  "special  solutions".  Such  solutions  do  involve 
t&e  complete  fluid  dynamic  equations  of  flow,  are  self- 
similar  in  nature,  and  require  that  the  initial  conditions 
be  precisely  those  which  insure  se’ f-similarity  (Ref  7:1). 

The  conversion  of  chemical  energy  of  the  propellant  to 
kinetic  energy  cf  the  shell  is  a  complex  process,  and  an 
exact  analytical  description  of  th?s  process  does  not  exist. 
Every  gun  internal  ballistics  analysis  employs  a  number  of 
simplifying  assumptions  which  reduce  the  problem  .o  a  model 
which  can  be  more  readily  analyzed.  In  order  to  provide  a 
basis  of  comparison  between  the  present  work  and  other 
studies  the  most  coamonly-used  assumptions  found  in  other 
gun  ballistics  studies  will  now  be  discussed. 

One  of  the  first  simplifications  applied  to  the  gun 
problem  was  that  the  propellant  was  completely  burned  before 
the  shell  motion  was  permitted  to  begin.  This  assumption 
was  used  by  Robins  and  Lagrange  in  the  eighteenth  century, 
then  by  Love  and  Pidduck  during  the  early  part  of  the 
twentieth  century  (Ref  3:347),  and  finally  in  a  more  modern 
work  by  Seigel  (Ref  7).  A  second  assumption  includes  the 
presence  of  burning  propellant  in  the  analysis,  but  requires 
that  the  propellant  velocity  be  the  same  as  the  gas  velocity. 
The  works  cf  Baer  and  Hitchcock  are  examples  of  studies 
which  employ  this  assumption  (Refs  1,  2,  and  S) .  A  third 
assumption  is  that  the  chamber  may  bo  represented  by  an 
"effective  chamber"  which  has  the  same  volume  as  the  actual 
chamber  but  a  diameter  equal  to  that  of  the  bore.  Seigel 
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states  that  the  "effective  chamber"  assumption  is  invalid, 
however,  in  that  an  analysis  incorporating  this  assumption 
neglects  cer-ain  significant  compression  effects  which  occur 
as  a  result  of  the  area  change  in  the  chamber  (Ref  7:28). 

A  final  assumption  is  that  the  gas  density  is  only  a  function 
of  time.  The  works  of  Heiney  (Ref  4:5)  and  Hitchcock 
(Ref  5:4)  illustrate  that  this  assumption  is  a  sufficient 
condition  fcr  a  linear  gas  velocity  distribution. 

Such  assumptions  as  the  ones  discussed  above  do  serve 
to  reduce  the  gun  problem  to  a  more  ammenable  form,  but  in 
doing  so  they  tend  to  fora  models  which  deviate  somewhat 
from  physical  reality.  One  might  well  question  the  validity 
of  a  model  which  represents  a  chambered,  propellant-burning 
gun  with  ^  non-chambered ,  non-propellant  burning  tube,  yet 
the  results  of  many  such  simplified  analyses  agree  quite 
satisfactorily  with  experimental  results  (Refs  1,  2,  4,  5, 

6  and  7).  Recently,  however,  evidence  has  appeared  which 
indicates  that  conventional  theory  is  not  always  providing 
acceptable  results,  particularly  for  power-burning  guns  with 
muzzle  velocities  in  excess  of  5000  feet  per  second  (fps). 
Baer  points  out  that  as  muzzle  velocities  '.icrease  beyond 
5000  fps  conventional  internal  ballistic  tucory  is  unable  to 
predict  detailed  gun  performance  (Ref  1:535).  Further, 
Hitchcock  noted  an  increased  deviation  between  theory  and 
experiment  for  muzzle  velocities  in  the  5500  -  6500  fps 
range  (Ref  5:25-26), 
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The  Present  Work 

The  present  work  has  two  objectives.  The  'irst  is  to 
model  the  gun  internal  ballistics  process  in  such  a  manner 
as  to  be  independent  of  the  four  commonly-used  assumptions 
described  above.  The  second  objective  is  to  use  the  results 
of  this  model  to  explain  the  deviations  between  theory  and 
experiment  for  high-speed  guns  noted  by  Baer  and  Hitchcock. 

Section  II  is  a  two-part  section  devoted  to  a  discussion 
of  the  analytical  model  used  in  this  work.  The  first  part 
contains  a  discussion  of  the  fundamental  assumptions  used  as 
a  basis  for  the  model,  while  the  second  part  contains  a  word 
description  of  the  operation  of  the  model.  Section  III 
contains  the  working  equations,  the  derivation  of  these 
equations,  and  a  detailed  list  of  the  assumptions  used  in 
deriving  them.  Section  IV  contains  a  discussion  of  the 
results.  The  conclusion  reached  as  a  result  of  this  work 
and  some  recommendations  for  future  work  are  contained  in 
Section  V, 

A  Fortran  Extended  computer  program  was  created  to 
perform  the  numerous  storing,  searching,  and  computational 
routines  involved  in  the  solution.  Pertinent  information 
about  the  type  of  computer  used,  storage  requirements, 
program  run  time,  as  well  as  a  simplified  logic  diagram,  a 
program  listing,  and  a  sample  output  are  contained  in  the 
appendices,  A  brief  glossary  of  terms  peculiar  to  the  gun 
internal  ballistics  field  is  proyided  in  Appendix  E. 
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1 1 .  The  Analytical  Model 


Fundamental  As  sumptions 

The  analytical  model  is  based  upon  two  fundamental 
assumptions.  The  first  is  that  for  small  but  finite  incre¬ 
ments  of  time  the  actual  internal  ballistics  process  may  be 
represented  by  an  artificial  sequence  of  four  separate 
"phases".  This  sequence  consists  of  a  constant  volume 
energy  transfer  phase,  a  shell  motion  and  wave  propagation 
phase,  a  propellant  motion  phase,  and  a  gas  expansion  and 
mass  transfer  phase.  The  implication  of  the  assumption  is 
that  for  small  but  finite  increments  of  time  the  net  result 
of  this  artificial  sequence  of  phases  is  approximately  the 
same  as  if  all  the  phases  had  occurred  simultaneously.  This 
sequence  bears  a  general  resemblance  to  a  thermodynamic 
cycle  in  that  a  system  is  changed  from  an  initial  set  of 
conditions  to  a  final  set  of  conditions  by  an  orderly 
progression  of  events;  for  this  reason  an  individual  sequence 
will  hereafter  be  referred  to  as  a  "cycle". 

The  second  fundamental  assumption  is  that  the  column  of 

t 

gas  and  propellant  between  the  breech  and  the  shell  may  be 
represented  by  a  fixed  number  of  individual  gas  "subvolumes" . 
At  any  instant  of  time  the  gas  properties  within  an  individual 
subvoluae  are  considered  to  be  constant.  These  property 
values  may,  however,  vary  from  one  subvolua.e  to  another. 

Gas  and  propellant  mass  transfer  may  occur  between  subvolumes, 
but  only  at  separate  and  specified  times  during  the  cycle. 
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The  purpose  of  these  two  assumptions  is  to  simplify  the 
internal  ballistics  process  into  one  that  is  more  readily 
analyzed.  The  first  assumption  separates  the  complex 
internal  ballistics  process  into  more  fundamental  processes: 
constant  • olume  combustion,  one-dimensional  motion  and  mass 
transfer,  and  finite-amplitude  wave  propagation.  The  second 
assumption  simnlifi»s  the  analysis  of  the  gas  and  propellant 
column  by  separating  it  into  a  number  of  small  constant- 
property-value  subvolumes.  These  subvolumes  are  then  analyzed 
using  the  theory  of  finite-amplitude  waves  in  gases,  a  basic 
energy  equation,  and  simple  equations  of  mass  motion  and  mass 
transfer. 

Word  Description  of  the  Model 

Overall  Physical  Description .  The  diameter  change  from 
the  chamber  to  the  bore  normally  occurs  over  a  finite  length 
of  the  gun  barrel.  For  this  analysis  this  are.',  change  is 
considered  to  occur  at  a  single  location.  No  other  changes 
in  the  physical  description  of  the  gun  barrel  are  made.  In 
this  work  the  region  between  the  breech  and  the  area-change 
location  will  be  referred  to  as  the  "chamber",  while  the 
remainder  of  the  gun  barrel  will  be  referred  to  as  the  "bore". 

The  projectile  is  initially  positioned  at  the  location 
of  the  area  change  and  the  chamber  is  divided  into  a  fixed 
number  of  cylindrical  segments.  The  axis  of  rotation  of 
each  segment  is  the  same  as  the  axis  of  rotation  of  the 
chamber.  All  segments  initially  contain  the  same  quantity 
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of  gas  mass,  the  same  number  of  propellant  grains,  and  have 
the  same  volume.  Further,  all  segments  initially  have 
identical  gss  property  values. 

The  Ga.<;  Subvolumes.  The  gas  which  occupies  the  avail¬ 
able  space  between  the  boundaries  of  a  single  segment  forms 
the  gas  sitbvolume.  Gas  subvolume  properties  change  as  a 
result  of  expansion  (wave  propagation)  and  mass  change. 

Gas  mass  change  occurs  as  a  result  of  mass  transfer  across 
the  gas  boundaries  and  as  a  result  of  a  propellant  burn. 

The  Propellant  Segments.  The  amount  of  mass  released 
by  a  given  mass  of  burning  propellant  during  a  finite 
increment  of  time  is  dependent  upon  the  surface  area  of  the 
propellant,  the  relative  flow  of  gas  past  the  propellant, 
and  the  pressure  of  the  gas  surrounding  the  propellant.  In 
order  to  account  for  the  surface  area  of  the  propellant  as 
it  burns,  the  number  of  propellant  grains  in  ^ach  propellant 
segment  is  fixed  at  the  initial  value.  Also,  all  grains 
within  a  single  propellant  segment  are  considered  to  burn  at 
the  same  rate.  Hence  all  the  grains  within  a  single  propel¬ 
lant  segment  are  kept  identical  with  one  another,  and  the 
mass  released  by  a  single  propellant  segment  during  a  single 
burn  time  increment  is  simply  the  mass  released  by  a  single 
grain  times  the  number  of  grains  in  the  segment.  Since  the 
number  of  grains  in  a  single  segment  is  fixed,  it  can  be  seen 
that  the  mass  of  propellant  in  a  single  segment  can  only 
decrease.  Location  of  the  various  propellant  segments  is 
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accomplished  by  fixing  the  length  of  each  segment  at  the 
original  value. 

There  is  no  requirement  that  a  propellant  segment  be 
located  entirely  within  a  gas  subvoluae;  therefore  if  a 
single  propellant  segment  happens  to  be  located  such  that 
its  length  is  divided  by  a  gas  boundary  the  burning  rate  of 
that  segment  should  actually  be  influenced  by  the  two  differ¬ 
ent  velocities  and  pressures  of  the  two  subvolumes.  When 
this  situation  occurs  the  burning  rate  of  the  entire  propel¬ 
lant  segment  is  determined  by  the  average  pressure  and 
relative  velocity  of  the  two  subvoluaes. 

The  Gas  Boundaries .  The  gas  boundaries  have  three 
functions.  The  first  is  to  serve  as  solid,  fixed  boundaries 
during  the  constant-volume  energy  transfer  phase  of  the 
cycle.  The  second  function  is  to  act  as  planes  of  mass 
transfer  during  the  gas  expansion  and  propellant-motion 
phases  of  the  cycle.  The  third  function  is  to  serve  as 
locations  for  the  finite-amplitude  waves  which  are  used  to 
change  the  gas  property  values  following  an  incremental  shell 
motion.  If  a  wave  travelling  toward  the  breech  is  designated 
an  "upstream’'  (against  the  jas  flow)  wave  and  a  wave  trav¬ 
elling  toward  the  shell  is  designated  a  "downstream"  wave, 
then  it  can  be  seen  chat  there  will  be  four  possible  types 
of  waves:  an  upstream  expansion  wave,  an  upstream  compression 
wave,  a  downstream  expansion  wave,  and  a  downstream  compres¬ 
sion  wave.  Gas  properties  are  changed  whenever  one  of  these 

a 

waves  travels  across  a  gas  subvolume.  Since  the  gas 
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properties  within  each  subvolurae  are  required  to  be  constant 
at  any  instant  of  time,  a  wave  is  not  permitted  to  be 
located  between  two  gas  boundaries.  A  wave  is  propagated  to 
the  next  boundary  only  if  it  is  determined  that  there  is 
sufficient  time  remaining  in  the  time  increment  for  this 
motion  to  occur,  if  the  wave  cannot  cross  at  least  half-way 
across  the  subvolume  it  is  fixed  at  its  current  location. 

Gas  Subvolume  Containing  the  Area  Change.  The  gas  sub» 
volume  which  contains  the  area  change  from  the  chamber  to 
the  bore  is  an  exception  to  the  normal  constant-diameter  gas 
subvolumes,  and  is  therefore  treated  somewhat  differently. 

A  "property  discontinuity*'  is  considered  to  exist  at  the 
location  of  the  area  change,  and  the  gas  property  values  on 
the  chamber  side  of  the  subvolume  are  not  necessarily  the 
same  as  the  gas  property  values  on  the  bore  side.  Further, 
unlike  other  subvoiumes ,  internal  gas  mass  transfer  does 
occur  from  the  chamber  side  to  the  bore  side.  This  process 
is  described  in  the  Word  Description  of  the  Model  Operation 
section  below. 

Word  Description  of  the  Model  Operation .  In  order  to 
illustrate  the  operation  of  the  model  as  well  as  clarify  the 
functions  of  the  previously-described  features  a  brief 
description  of  a  typical  cycle  will  now  be  given. 

The  first  phase  of  the  cycle  is  the  constant-volume 

energy  transfer  phase.  All  motion  is  frozen,  then  each 

propellant  segment  is  burned  for  a  single  time  increment. 

«  ♦ 

This  burning  process  changes  the  pressure,  density. 
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temperature,  and  gas  volume  in  each  subvolume.  The  mass  of 
each  propellant  segment  is,  of  course,  reduced.  After  this 
increment  of  burning  is  completed  and  the  appropriate  adjust¬ 
ments  have  been  made  in  the  subvolume  property  values  the 
cycle  proceeds  to  the  second  phase. 

The  second  phase  is  the  projectile  motion  and  wave 
propagation  phase.  The  average  of  the  pressure  before  the 
burn  and  the  pressure  after  the  burn  in  the  subvolume  adjacent, 
to  the  projectile  is  considered  to  act  upon  the  shell  for  one 
time  increment.  This  force  changes  the  projectile  accelera¬ 
tion  and  velocity  and  moves  the  shell  to  a  new  location. 

The  change  in  projectile  velocity  is  considered  to  be  impul¬ 
sive  and  to  generate  a  finite  expansion  wave  which,  later 
in  the  cycle,  will  propagate  toward  the  breech.  The  model 
first  propagates  expansion  waves  from  earlier  projectile 
motion,  then  procedes  *o  other  types  of  waves  which  may  be 
present,  including  the  recently-generated  expansion  wave 
from  the  current  projectile  motion.  If  a  compression  wave 
"catches  up"  with  another  wave  of  like  kind  and  direction 
the  wave  strengths  are  combined  prior  to  further  propagation. 
Different  types  of  waves  or  like  waves  travelling  5n  opposite 
directions  are  not  combined.  Waves  which  encounter  the 
projectile  or  breoch  are  reflected  in  like  kind,  and  s  wave 
which  crosses  the  gas  subvolume  containing  the  area  change  is 
reflected  as  two  waves  of  appropriate  strength  and  direction 
(Ref  7:28).  When  all  present  waves  have  been  propagated  as 


10 


GAM/ME/72-2 


far  as  possible  during  a  single  time  increment  the  cycle 
proceeds  to  the  third  phase. 

The  third  phase  of  tho  cycle  is  the  propellant  motion 
phase.  The  pressure  drop  across  the  length  of  the  segment 
is  determined,  then  an  estimated  drag  coefficient  and  an 
estimated  effective  area  are  used  in  r.  simplified  equation  of 
motion  to  determine  the  new  velocity  and  position  of  the 
propellant  segment.  When  all  the  segments  have  been  moved 
the  cycle  proceeds  to  the  fourth  and  final  phase. 

The  last,  phase  of  the  cycle  is  the  gas  expansion  and 
mass  transfer  phase.  The  net  effect  of  each  wave  that  has 
crossed  a  ingle  gas  boundary  is  used  to  determine  the  new 
velocity  of  that  boundary.  Once  the  new  velocity  is  deter¬ 
mined  the  boundary  is  moved  at  that  velocity  for  a  single 
time  increment.  The  boundary  velocity,  barrel  cross- 
sectional  area,  and  the  gas  density  of  the  next  downstream 
subvolume  are  used  to  determine  the  gas  mass  transfer  across 
the  boundary  during  this  motion.  After  all  boundaries  have 
been  relocated  the  total  propellant  ma^s  ,/ithin  each  sub¬ 
volume  is  redetermined.  Kith  the  su’/volume  pressure  held  at 
the  value  determined  during  the  wave  propagation  phase,  the 
remainder  of  the  gas  property  values  arc  then  determined  from 
an  equation  of  state.  This  final  property  determination 
marks  the  end  of  the  cycle. 

The  cycle  just  described  is  repeated  until  it  is  deter¬ 
mined  that  the  projectile  pos.  ion  exceeds  the  length  of  the 
barrel,  at  which  time  the  analysis  ends.  Provisions  are  made 
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in  the  Model  to  check  for  propellant  burn-out  prior  to 
projectile  exit  in  order  that  the  propellant-burning  and 
motion  parts  of  the  cycle  may  be  deleted.  Figure  1  illus 
tr&tes  the  physical  appearance  of  the  model  prior  to 
projectile  motion  and  at  some  later  time. 
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III.  Analysis 


Analytical  Assumptions 

The  two  fundamental  assumptions  which  form  the  basis  for 
the  analytical  model  have  been  previously  discussed  in  the 
first  part  of  Section  II.  The  following  analytical  assump¬ 
tions  have  also  been  used; 

(a)  All  motion  is  one-dimensional. 

(b)  Propellant  burning  takes  place  under  constant- 
volume  conditions. 

(c)  The  propellant  burning  rate  is  a  function  of  gas 
pressure  and  relative  gas-to-propellant  velocity. 

(d)  The  propellant  grains  burn  uniformly  over  their 
entire  surface. 

(e)  The  gas  obeys  the  Nobel-Abel  equation  of  state  with 
a  constant  covolume. 

(f)  The  gas  has  a  constant  ratio  of  specific  heats. 

1 

(g)  The  gas  subvolume  boundaries  are  adiabatic. 

(h)  The  drag  coefficient  for  the  propellant  segments  is 
constant . 

(i)  The  drag  on  the  projectile  consists  of  a  variable 
aerodynamic  drag  and  a  constant  friction  drag. 

(j)  The  diameter  change  from  the  chamber  to  the  bore 
occurs  at  a  single  location. 

(k)  A  normal  shock  forms  ahead  of  the  projectile  as 
soon  as  the  projectile  metier-  begins. 
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(1)  The  presence  of  the  propellant  exerts  no  influence 
on  the  wave  propagation  process. 

The  following  are  considered  to  be  negligible: 

(a)  Heat  transfer  to  the  gun  walls  and  to  the  projectile. 

(b)  Friction  losses  between  the  gas  and  the  gun  walls. 

(c)  Friction  losses  between  the  propellant  and  the  gun 
walls . 

(d)  Drag  due  to  projectile  rotation  (nfling  drag). 

(e)  Losses  due  to  propellant  gas  leaking  past  the 
projectile . 

(f)  The  pressure  gradient  between  the  front  of  the 
projectile  and  the  downstream  side  of  the  norisal 
shock. 

i  (g)  Effects  due  to  gun  recoil. 

(h)  Effects  due  to  variations  in  the  initial  teeperature 
of  the  propellant. 

The  Working  Equations 

The  Energy  Equation.  The  first  law  cf  thermodynamics 
for  a  constant  volume  subvolume  with  no  mass  flow  is 

Q  ■  AU  (1) 

where,  for  the  constant  volume  adiabatic  combustion  used  in 
this  work 

q  .  — L-  (Ref  3:175)  (2) 

TS 

and 
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Au  *  CyAT 


-here  Yg  ,  propellant  gas  ratio  of  specific  heat: 
For  the  gun  problem 


Cv  E  V — rr  (Ref  8:126) 

fi 


Since  P,  the  “force  constant’*  is  defined 


F  =  R  T. 

g  iso 


Eq  (1)  becomes 


Vise  .  V«) 

Y*  -  1  Yg  '  1  (6 

For  a  finite  quantity  of  Eq,  (2)  and  (S)  „y  k.  expressed 


R  T 

Q  a  AM 

YS  '  1  * 


AU  s 


•H  'Vf  -  Vi> 


where  *  final  mass  of  gas 

Mi  *  initial  mass  of  gas 
«  final  gas  temperature 
Ti  *  initial  gas  temperature 

Equating  (7)  and  (8) 


S  T 

— j B — liS.  fu  m  u  \ 

yg  -  i  V 
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T 


f 


S 


T. 

ISO 


C  i 


♦  T. 
1 


M. 

x 


(10) 


Equation  of  State .  The  equation  of  state  used  is  the 
"Nobcl-Abcl"  equation  of  state  with  a  constant  covolume 

P (V  -  Mb)  a  MR  T  (11) 

where  b  «=  covolune 


Mass  Change  Due  to  Propellant  Burn .  The  change  in  mass 
for  a  single  propellant  segment  during  a  single  time  increment 
is 


where 


M 

P 

il 

g 


change  in  a  single  grain  volume 

propellant  density 

number'  of  grains  per  segment 

propellant  mass 

gas  mass 


(12) 


The  propellant  burn  rate  is  taken  to  be 


R  =  3(P/10C0)n  +  K  v  (Ref  5:9)  (15) 

where  R  =  propellant  burn  rate  (length/time) 

3  =  burn  rate  at  1000  psia  and  v?  c  0 
P  3  gas  pressure 
n  =  burn  rate  exponent 

K  =  erosive  burn  constant 
c 

v?  =  relative  gas-to-propcl lant  velocity 
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The  volume  change  of  a  single  grain  is 


4Vs "  "W  '  W 


where  Ap  «  grain  surface  area 
Rp  »  grain  radius 

(  )’  *  value  after  proo*li^nt  burn 


The  absolute  value  in  Fq  (14)  is  necessary  because  some 
grains  are  designed  such  that  the  surface  area  increases 
during  the  initial  burn  process.  The  surface  area  vs.  grain 
racius  for  the  particular  problem  studied  was  obtained  from 

t 

tabular  data  (Ref  5:38-41).  Equations  (12),  (13),  and  (14) 
are  used  to  determine  vs.  M.  in  Eq  (10). 


Wave 


:ion.  The  pressure  change  induced  by  a 


finite-amplitude  wave  of  strength  Atf  is  given  by 


AP  -  -pa(Av)  (Ref  7:10-12)  (IS) 


where  AP  *  finite  pressure  change 

p  *  gas  density  ahead  of  wave 
a  «*  sonic  velocity  ahead  of  wave 
Av  ■  finite  velocity  change 


Wave  velocity  is  given  by 


a  t  v  (Ref  7 ; 1 1 ) 
8 


(16) 


where  v^  »  wave  velocity 
v  »  gas  velocity 
a  «  gas  sonic  velocity 
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As  mentioned  previously  in  Section  II,  a  wave  which  encounters 
the  gas  subvolune  containing  the  chamber-bore  area  change  is 
split  into  two  waves.  For  example,  a  bore-side  upstream¬ 
travelling  expansion  wave  is  split  into  a  compression  wave 
which  travels  back  toward  the  shell  and  an  expansion  wave 
which  continues  on  toward  the  breech.  The  strengths  and 
directions  of  the  split  waves  are  determined  in  the  following 
aauner . 

Consider  an  upstream-travelling  rarefraction  (expansion) 
wave  AVj  which  has  just  reached  the  bore  side  of  the  gas 
subvolume  containing  the  chamber-bore  area  change  (Fig.  2a). 
The  change  in  pressure  on  the  bore  side  of  the  subvclume  is 
determined  with  Eq  (IS) 

-  -PbVvJ 

where  pb  *  bore-side  density 

afa  «•  bore-side  sonic  velocity 

The  wave  is  advanced  to  the  point  of  area  change  and  the  new 
bore-side  pressure  and  velocity  values  are  determined 
(Fig.  2b) 


%  *  Api 


(17) 


Vb  *  Arl 


(18) 


The  decreased  pressure  on  the  bore  side  induces  as  increased 
mass  flow  from  the  chamber  side.  The  amount  of  mass 


Id 
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transferred  fro*  the  chamber  side  to  the  bore  side  is 


'  0chVb1(4t> 

where  pr^  ■  chamber  side  density 
■  bore  side  area 
At  *  time  increment 


(19) 


The  amount  of  gas  aass  on  the  bore  side  is  therefore 
increased  by  an  amount  AM^  while  the  gas  mass  on  the  chamber 
side  is  decreased  by  the  same  aaount  (Fig.  2c).  If  tempera¬ 
ture  is  assumed  to  be  constant  during  this  process  the  new 
bore-side  pressure  becomes 


(«b  ♦  MPV,, 
\  ♦  -Zh.H 

while  the  new  chaaber-side  pressure  is 

(M.  -  AM  )R  T 

P  .  _J_b  _ .Ll-g  . 

Cj  (Vc  -  b(Mc  -  AMj)] 


where  *  bore-side  teaperature 

Tc  *  chamber-side  temperature 


(20) 


(21) 


It  can  be  seen  from  Eqs  (20)  and  (21)  that  will  be 

greater  than  P.  (but  still  less  than  P.)  while  P  will  be 
bl  b  C1 

less  than  Pc»  Hence  th*  net  effect  is  to  produce  a  downstream- 

travelling  compression  wave  of  strength  AP«  «  P^  - 
(positive)  at  Aa  (Fig.  2a)  and  an  upstream-travelling  expan¬ 
sion  wave  of  strength  AP,  »  P  *  P  (negative)  et  B  in 

a  c,  c 
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in  Fig.  Id.  A  similar  analysis  holds  tor  the  other  types 
of  waves. 

Propellant  Segment  Motion .  Propellant  segment  motion 
is  determined  from  a  simplified  equation  of  motion.  It  is 
assumed  that  the  sum  of  the  forces  Ef^  acting  on  a  single 
propellant  segment  for  a  single  time  increment  is 

Ef  -  (AP)  A  ♦  D  (22) 

pep 

where  AP  *  pressure  difference  across  the  segment  length 

A  *  estimated  "effective  area" 
e 

Dp  *  estimated  aerodynamic  drag 
The  estimated  effective  area  of  the  segment  is  taken  to  be 

V 

A  •  r2-  (23) 

P 


where  *  volume  of  the  propellant  segment 

Lp  ®  fixed  length  of  the  propellant  segment 


The  estimated  aerodynamic  drag  on  the  propellant  segment  is 
taken  to  be 


D  *  1/2  p  v  A  C. 
p  g  r  e  d 


(24) 


where  *  gas  density 

v„  ■  relative  velocity  of  gas  past  tho  propellant 
•  estimated  drag  coefficient  (constant) 


If  the  acceleration  of  the  segment  a^  is  approximated  by 
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where  P^  *  aerodynamic  drag  pressure 

P  *  ambient  (upstream  of  shock)  pressure 

u 

ym  «  ambient  ratio  of  specific  heats 
v  *  projectile  velocity 
a  *  ambient  sonic  velocity 


Projectile  Equation  of  Motion,  The  equation  of  motion 
fpr  the  projectile  is  taken  to  be 


p,d>At 


)  c„) 

pr 


(30) 


where 


Ppr  =  pressure  on  the  breech  side  of  the  projectile 
Pf  «  estimated  constant  friction  pressure 
Pa<l  *  aerodynamic  drag  pressure 
M  *  projectile  mass 
a pr  «  projectile  acceleration 


24 


mamai — - — x&maaSsiasJ&ai 


sum*— 


GAM/ME/72-2 

IV.  Discussion  of  Results 

Of  the  two  objectives  stated  in  the  introduction  to 

this  work  only  the  first  has  been  met,  and  that  with  extremely 

limited  results.  The  model  has  performed  acceptably  for  only 

one  complete  cycle  following  the  initial  motion  of  the  shell. 

Results  obtained  from  the  second  cycle  indicate  that  the  model 

is  failing  to  combine  gas  and  propellant  motion  in  such  a 

manner  as  to  obtain  a  realistic  pressure  distribution. 

Specifically,  an  unrealistically  low  pressure  in  the  gas 

subvolume  adjacent  to  the  shell  is  established  early  during 

the  second  cycle.  This  low  pressure  in  turn  sets  up  wave 

reflections  from  the  chamber  subvolume  with  unrealistically 

high  velocity  strength  values.  These  erroneous  waves  are 

then  propagated  for  the  remainder  of  the  time  increment, 

and  the  resulting  model  bears  little  resemblance  to  the 

actual  physical  situation.  The  excessively  low  pressure  is 
* 

a  direct  result  of  the  separation  of  gas  motion  from  propel¬ 
lant  motion.  The  sudden  expansion  of  the  subvolume  adjacent 
to  the  shell  caused  by  the  shell  motion  increases  the  volume 
available  to  the  gas  in  that  subvolume.  Because  there  is  no 
motion  of  propellant  into  that  subvolume  at  that  point  in 
the  cycle,  the  space  that  should  be  occupied  by  some  portion 
of  propellant  mass  is  not;  hence  the  gas  expands  to  fill  an 
unrealistically  large  volume  and  the  pressure  drops. 
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V.  Conclusion  and  Recommendations 


Conclusion 

The  conclusion  reached  as  a  result  of  the  work  to  date 
is : 

The  basic  model  concept  shows  promise,  but  at 
present  the  model  is  failing  to  realistically  represent  the 
combined  flow  of  gas  and  propellant. 

Recommendations 

The  following  recommendations  are  made  concerning  this 
analysis : 

1.  The  basic  model  concept  should  be  revised  to 
realistically  represent  the  combined  flow  of  gas  and  propel¬ 
lant  in  order  to  preserve  a  realistic  pressure  gradient. 

2.  The  propellant  segment  drag  coefficient  should  be 
computed  from  Reynold's  number  considerations  instead  of 
the  present  constant  value. 

3.  A  more  accurate  representation  of  the  shell  sliding 
friction  should  be  attempted. 

4.  An  allowance  for  energy  loss  due  to  heat  transfer 
should  be  introduced. 

5.  Propellant  burning  rate  values  should  be  obtained 
from  tabular  data  rather  than  the  conventional  pressure  and 
erosive  burn  scheme  used  in  this  work. 


26 


GAM/ME/72-2 


Bibliography 


1.  Baer,  ?.  G.  "Prediction  of  High  Velocity  Solid  Propellant 
Gun  Performance  by  Gas  Dynamic  Computer  Program." 

Ballistic  Research  Laboratories  Report,  Aberdeen  Proving 
Ground,  Maryland.  1970. 

2.  Baer,  P.  G„  "A  Mass  Point  Computer  Program  for  the  Gas 
Dynamic  Problems  of  High  Velocity  Interior  Ballistics." 
Unpublished  report  by  the  Ballistic  Research  Laboratories, 
Aberdeen  Proving  Ground,  Maryland.  Undated. 

3.  Corner,  J.  Theory  of  the  Interior  Ballistics  of  Guns . 

New  York:  Johr  Wiley  and  Sons,  Inc.  1950. 

4.  Heiney,  D.  K.  "Simplified  Interior  Ballistics  of  Closed 
Breech  Guns."  AFATL-TR-67-42 .  Air  Force  Armament 
Laboratory  Report,  Eglin  Air  Force  Base,  Florida. 

April  1967. 

5.  Hitchcock,  J.  E.  "Gun  Interior  Ballistics."  AFRPL-TR-69-21 1 . 
Air  Force  Rocket  Propulsion  Laboratory  Report,  Edwards  Air 
Force  Base,  California.  September  1969. 

6.  Seigel,  A.  E.  "The  Theory  of  High  Speed  Guns."  NATO 
Agardograph  91.  North  Atlantic  Treaty  Organization 
Advisory  Group  for  Aerospace  Research  and  Development, 

United  States  Naval  Ordnance  Laboratory,  White  Oak, 

Silver  Spring,  Maryland.  May  1965. 

7.  Thornhill,  C.  K.  "A  New  Special  Solution  to  the  Complete 
Problem  of  the  Internal  Ballistics  of  Guns."  RARDE  Report 
2167.  Royal  Armament  Research  and  Development  Establish¬ 
ment,  Fort  Halstead,  Kent,  England.  February  1967. 


27 


taSei 


saxaul 


—Sag! 


GAM/KE/72-2 


Appendix  A 

Computer  Program  Features  and  Requirements 
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Computer  Program  Features 

Debug  Mode  of  Operation .  A  self-debugging  feature  is 
built  into  the  program  to  enable  automatic  debugging  during 
future  program  modifications.  The  self-debugging  feature 
i’s  activated  by  replacing  the  "FTN.1'  control  card  at  the 
beginning  of  the  deck  with  an  ”FTN(D)“  control  card,  and 
increasing  the  memory  requirement  by  12K.  The  debug  feature 
causes  the  following  to  occur  automatically: 

a.  Automatic  bounds  check  on  all  arrays. 

b.  Printout  of  certain  key  program  variables 
values  whenever  these  value,  change,  along 
with  the  program  location  of  the  change. 

Logic  tracing  is  available  by  adding  a 

C$  TRACE 

col:  12  7 

card  immediately  after  the  "C$  DEBUG”  card  in  the  deck. 
Further  information  on  the  debug  mode  of  operation  is  con¬ 
tained  in  Chapter  11  of  the  Control  Data  6400/6500/6600 
Computer  Systems  Fortran  Extended  Reference  Manual. 

Solution  of  the  Preburned  Propellant  Problem.  The  pre¬ 
burned  propellant  problem  may  be  considered  with  this 
program  by: 

a.  Entering  all  propellant  data  as  if  the 
propellant  were  going  to  be  burned. 

b.  Setting  the  value  of  the  'variable  ‘’NOCH"  to 
"2"  on  the  appropriate  dsta  card. 
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Computer  Program  Requirements . 


Language : 
Computer: 
Storage 


Run  time: 


Fortran  Extended. 

Control  Data  6600  (Digital) 
36K  (Binary) 

4SK  (Compile/no  debug) 

60K  (Debug) 

Undetermined. 


GAM/MS/72-2 


Comput 

NOTE: 

Symbol 

A 

ADP 

AGRO 

AGV 

A1 

BA 

BD 

BETA 

BGMAS 

BP 

BRO 

BT 

BURNA 

BV 

BXP 

CA 

CD 

CGMAS 

CHD 

CHL 

CHRO 


er  Program  Symbols 

The  term  "bore- side"  refors  to  the  bore  side  of  the 
subvolumc  containing  the  chamber-bore  area  change. 
"Chamber-side"  refers  to  the  cnamber  side  of  the  same 
subvolume. 

Definition 

Area 

Aerodynamic  drag  pressure 
Average  gas  density 
Average  gas  velocity 
Multipurpose  variable 
Bore  area 

Bore  diameter  (input) 

Pressure  burn  coefficient  (input) 

Bore-side  gas  mass 

Bore-side  pressure 

Bore-side  gas  density 

Bore-side  temperature 

Propellant  burn  area  (tabular  input) 

Bore-side  gas  velocity 
Pressure  burn  exponent  (input) 

Chamber  area 

Propellant  drag  coefficient  (input) 
Chamber-side  gas  mass 
Chamber  diameter  (input) 

Chamber  length  (input) 

•  Chamber-side  gas  density 


32 


GAM/ME/72“2 


Symbol 

CHV 

CL 

CMAS 

CMASS 

CM  I G 

COVOL 

CP 

CRO 

CT 

CV 

CX 

DCM 

DIST 

DM 

DMC 

DHL 

DMR 

D? 

DT 

DV 

DVB 

DVC 

DVSUM 

EBK 

F 


Definition 

Chamber-side  gas  velocity 

Fixed  length  of  propellant  segments 

Propellant  segment  mass 

Total  propellant  mass  in  subvolume  (output  only) 

Igniter  charge  mass  (input) 

Covoluae  (input) 

Chamber-side  pressure 

Propellant  density  (input) 

Chamber-side  temperature 

Propellant  segment  velocity 

Propellant  segment  position 

Mass  Change  due  to  burn 

Distance  (various  uses) 

Incremental  mass  change  (various  uses) 

Center  section  of  DCM  for  KTYP  «  10  and 
KTYP  -  11  (Ref  to  Fig.  3). 

Left  side  of  DCM  (Refer  to  Fig.  3). 

Right  side  of  DCM  (Rofer  to  Fig.  3). 

Pressure  change 

Time  increment 

Velocity  change 

Bore-side  velocity  change 

Chamber-side  velocity  change 

Sum  of  velocity  changes  at  an  individual  gas 
boundary 

Erosive  burn  constant  (input) 

Cun  constant  (also  known  as  "force  constant") 
(input) 
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JB 

-JC 

JS 

K 


Definition 

Friction  pressure  (input) 

Ambient  ratio  of  specific  heats  (input) 
Propellant  gas  ratio  of  specific  heats  (input) 
Subvolume  gas  mass 

Initial  propellant  grain  mass  (input) 

Number  of  grains  per  propellant  segment 
Grain  burn  radius  (tabular  input) 

Gun  barrel  length  (input) 

Counter  (various  uses) 

Bore  gas  boundary  reference 
Program  section  identifier 
Stored  gas  boundary  value 

•'Is  wave  available  for  propagation”  indicator 

IWA  ■  0  Wave  present  and  ready  for 
propagation 

IWA  "1  No  wave  present 

IWA  ®  2  Wave  present  but  already 

propagated  during  this  time 
increment 

Index  (used  during  wave  propagation) 

Wave  type;  J  *  1  Upstream  expansion 

J  «  2  Downstream  expansion 
J  «  3  Upstream  compression 
J  »  4  Downstream  compression 

Bore-side  wave  type 

Chamber-side  wave  type 

Stored  wave  type  value 

Counter  (various  uses) 
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Symbol 

KCHN 


KTYP 

L 

I 

LTYP 


M 

N 

NB 

NOCH 


NTAB 


OBP 

OCP 

OP 


ORAD 

OSHV 

OSURFA 

P 

PA 


Definition 


Index: 


KCHW 


0  Wave  is  not  chamber 
reflection 

KCHW  »  1  Wave  is  chamber 

reflection  from  area 
change  subvolume 

KCHW  »  2  Wave  is  bore  reflection 
from  area  change 
subvolume 


Type  of  propellant  segment  (Refer  to  Fig.  3) 
Counter  (various  uses) 


Index: 


LTYP  *  0  Propellant  segment 
within  subvolume 
LTYP  *  1  Propellant  segment 
divided  by  upstream 
gas  boundary 


Counter  (various  uses) 

Propellant  segment  counter 
Number  of  g#s  boundaries  (input) 


Index : 


NOCH  *  0  Propellant  segment 
present  (input)  ins 
subvolume 

NOCH  »  1  Propellant  segment  not 
present  in  subvolume 
NOCH  3  2  No  propellant  in  gur. 


Number  of  tabular  entries  in  the  grain  surfr.ee 
area  (input)  vs  burn  distance  table 


Bore-side  pressure  before  incremental  burn 
Chamber-side  pressure  before  incremental  burn 


Normal  subvolume  pressure  before  incremental 
burn 


Grain  uurn  radius  before  incremental  burn 
Shell  velocity  at  beginning  of  time  increment 
Grain  surface  area  before  incremental  burn 
Subvolume  gas  pressure 
Ambient  pressure  (input) 
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Symbol 

t’cmasb 

TCMASC 

TEMP 

TIME 

TISO 

TOTCM 

TTC 

TYPGUN 

TYPROP 

V 

YEL 

VOL 

W 

KDIS7 

WV 

X 

XMR 

XNB 

XRAO 


Definition 

Total  bore-side  propellant  mass 
Total  chamber-side  propellant  mass 
Temperature  (used  in  output) 

Time  expired 

Isochoric  flame  temperature  (input) 

Total  propellant  mass  (input) 

Time  to  cross  subvolume  (used  in  wave 
propagation) 

Type  of  gun  being  analyzed  (input) 

Type  of  propellant  (input) 

Subvolume  gas  velocity 
Velocity  (used  in  output) 

Volume  (various  uses) 

Stored  value  of  wave  strength 
Distance  (various  uses) 

Wave  velocity 

Gas  boundary  position 

Propellant  ma3s  ratio  (before  burn  vs  after 
burn) 

Same  as  NB 
Grain  burn  radius 
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Appendix  C 

Simplified  Computer  Program 
Logie  Diagram 
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Simplified  Cojjmter  Program  Logic  Diagram 

INITIALIZATION  SECTION 


GAM/ME/72-2 


SHELL  MOTION  SECTION 


* 
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WAVE  PROPAGATION  SECTION 


GAM/ME/72-2 


CHAMBER  WAVE  REFLECTION  SECTION 


DETERMINE 
STRENGTHS 
OF  EACH 
WAVE 


GAM/ME/72-2 


PROPELLANT  MOTION  SECTION 


44 


GAM/ME/72-2 


PROPELLANT  SEGMENT  LOCATOR 


SECTION 


SET  VALUE 
OF  VARIABLE 
*'NOCH"  TO  1 
(PROPELLANT 
NOT  PRESENT) 


Vypes  of  Propellant  Segment  Locations. 


GAM/ME/72-2 


Appendix  D 


Computer  Program  Listing 
and 

Sample  Owtpmt 
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GAM/ME/72-2 


CDC  6600  FTN  V 2.C-279C  OFT=^  2Z/J//72 


PROGPAM  GUN (INPUT, 0UT°UT, C E HUG = OUT PUT) 

C  Of  PUG 
CC  ap~AvP 

r  •'  ?TOPFS(nGVa,RF,P7,PV,CGVOl,CHV,CP,C7,CV,OIST  ,0F, 
CT  ST05Fc ( IWA ,IX,J,JS, JC,JF,K,KCHW,LtM,N,NOCH,P, FDIF 
Cr  STrOOr^(CHAK,OMA><,CVr,  nVO»CMK,ARATP,ARATT) 

C~  STORCS(T»TAv,TCMAS,TCMASB,TCMASC»TTC#V»VOL*W,X,CX 
C«  <sT0RrS(LTYP,KTYP,wniST,0VSUM,GMAS»8GMAStC5MAS,RC» 
C'  STOREF(OMR  ,CMAS> 


CV,I,  IE,  ID,  I,J 
,S  FA  ,SH V  ,SHX ) 


CHF0,6PC,0f',CHC) 


****■♦*****#*»**»♦ 


GUN  INTERNAL  BALLISTICS  PROBLEM 


OT^EnsICN  F(^0) ,V(FG) ,RO(50) ,T(50),GMAS(50),TCMAS 
OTHPNSION  CV  (50 ) ,CMAS (50) ,R<50) ,CXC50) 

OT MENTION  X  <S 1) ,  IWA(5t,4> ,W(5i,4) , BURKA (2 0) , G fi A D ( 


(5 f!)  ,DVSLM  (50 
EG) 

»*  *»*#*»***♦*♦*#* 


*  PFAU  INPUT  VALUES  * 

*  * 

orA0  POuO,TYPGUN,TYFRnP 
°  r  A  0  2G1i,GUNL,CHL,CH0,°D 

^tao  ?o°i,shm,totch,cmig 

RPAO  2  C  7  I , PSHCT , p  F 
PPAO  2J7i,F,TISO,GAflG,COVOL 
RfAO  P.71,CRO,GNHAS,CO 
READ  ?0'11»B2TA,RXP,FT< 

°r AO  ?J01,PA,TA,ROA,G5MA 
"‘AO  2^2,NTAP,NP 

?fi-D  2*03,  (GRAD(I)  ,BURNA(I)  ,i=i,ntapi 

=>rA0  2 1  j 2, NOCH 
o'HnTrl.O  JO. 

DT=. J  CG  7  5 

S V A=PGRT (1.4*32. 1 74*53. 35»TA) 
nT=.GGOCl 

♦  * 

♦  PPINT  INFUT  VAuUES  * 


Drpjx  i:jv- 
ppINT  3 j 10  »TYFGUN 
PMK1’  3j2r  , GUNL 
ppJNT  7J3:,CHL 
p r I N T  3ci'7,GwO 
PPINT  3-50, PQ 
ppINT  30En 
=PTNT  7  ;  ?'*  ,  p  >HOT 
Pf  INT  3  J  ^  j , FP 
P°I‘»T  3.PGfSHM 
PFINT  31  J'j 
PPIMT  3 1 10  , TYPROP 
ppI.N;  3120  ,T0TCN 


»***#* *********** 
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gun 


COC  6600  PTN  V2.0-279C  0FT=1  02/07/72  1 


3  P  T  M  t  313. *  CM IG 
3r  I  .‘I  T  3lL' ,r30 
PRINT  '*  1 6  i  t  T  I S 0 

PFTMT  ■?;«,'  ,F 

P f'  I N  T  3 1  7  ,  3  r  t  A 
33:nT  o„.'U1,nXP 

Dc>lKiT  3lO<"  ,  ff'K 

PRINT  ?2jT,C0VCL 
oc  I  NT  .?2i;,GAH4 
P'INT  3?2')  ,GNMAS 
P?IMT  32  3'  » CO 
PrIk,T  3?4r 
PRINT  3?r>:»PA 
33INT  TA 

PrlMT  •»??■'  ,  P.G  A 
PRINT  32P.:,SVA 
P0INT  3291 
PRINT  3 7  3  r  »  C T 
°PINT  3-»13,N« 

PRINT  7 3At 

PRINT  335r 

PRINT  336?  ,  (GPAO  (I)  ,9URNA  (I)  ,I  =  1,MAB> 


i  ■*  #** 


* 


* 


*  INITIALIZE  VARICJS  VA3IA3LFS  * 

*  LOG AT r  INITIAL  GAS  ROUNDARIES  ANO  CHANGE  SEGMENTS  * 

*  i  * 


PTr.  C-M'l 
xnp-  *;n 

T  T  M  c  =  p  T 

3r  =  p/iisc 

S  V  K= 0.  A  MG  *  3  2  •  1 74  *  RG 
3 A  =  .  TP^L'L'C*30 
CA  =  .  V#«54*CHC*CHD 
G^X~Cnl 

SHVrf-  . 

«uA=r. . 

CL=rwt  /(XNS-1.) 

GNS=TOTCM/ { (XNB-1.) *GNMAS> 

L-NP-2 

HIST-:  . 

PC  1 L  1  =  1,  L 
v  ( I )  =  R I S  T 
Ry(I)=CIST 
pvsiiN  (i )  =; . 

TIST~nI 3T+CL 
X  t.'.3)  =C^L 
nu  ?L  1=1 ,N3 
PO  ?u  J=i,4 
IWA (I, J) =1 
M  <  I  ,  J  )  =  “  .  • 


♦  .ITTCPNIUE  CONCITICNS  AT  PSK)T 


SO 


T- 1 
N  =  1 

C*as  ci >  =  TnTnM/c ynp-i, ) 

VOl=l?.*CA*X  C 2 ) -CMA^(i)/CRO 

GMASCt  >=V0t*R0A4-CKIG/  (XMR-i.) 

P(1)=32.*GMAS(1) *SG*TISO/{VOL-COVCL*GHAS(i>> 

T<1)=TIS0 

R(1)=GRAD(1) 

TOMAS  (l)=CMAS(i) 

AGV=0, 

A  C  V  =  0  • 

GO  TO  Z^t) 

30  IP(P(1)  .GF.PSHOT)  GO  TO  40 

TTHr=TIMt<CT 
GO  TO  200 
40  L=Mn-i 

n'  5  C  1-1,1 
(7  )  -  °  ( 1 ) 

VCI) =l  . 

PO(T)=RO(l) 

T(I)=TTSO 
G  VAS (I) =G«AS<1) 

TOMASCD  =TCMAS(i) 

ON  AS  (  T )  =CH  AS  (1) 
r-V(I)=0. 

50  °{I)=9(i) 

OP=o (i) 

T  =  *15—1 
OT=.LCd15 

SHA=22,i7U^(P(N0-l)-FP)*BA/SHK 
<U'V  =  S*A*OT 

SHX  =  S^X  +-SH V¥SHV/  ( SHA*2  • ) 

X  (Mn }  =SHX 

CT=TTSO 

OT=TTSO 

nv=SHV 

OV=5HV 

nr=-. 374 »PC( II *SCKT(SVK*TISO> *OV 
*>P  =  P  Cl)  +np 

t’GHAS=R0(I)*9A*3V*PT*tc« 

VOL- (SHX-CHL) ♦OA*i2e 
npprPr.Mfl^/vOL 
0“=PG’4  AS 


VOL=(CHL-X (I) }*CA*12.-CMASCD /CRO 
CGMAS=GM AS  (l)-OM 
CHRO=CGMAS/VOL 

0P=- • 374  *OHRO*SGR  T (SVK*TISO) *OV 
CP=P  (I>*QP. 

H  Cl, 1) =-2.68*0°/ (PO(l) *SORT(SVK*TISO) » 

IKA(I,1)=0 

CHV=W(I,1) 


51 


GAM/ME/ 72 -2 


CCC  6600  FTN  V3.C-270C  CFT  =  1  02/07/72  i 


7 CwA*?r=f'MA  S  (  i ) 

T  r  M  A  c> n  = ,'  t 
3VSUM  < I ) =W ( T  » 1 ) 
GO  TO  <*„’•> 


PROrfLL ANT  OURN  5E07ICN 


DTSTsX (Nn)-X(Nn-l) 

IF(X(7.-n-i)  .LT. CHL.ANO.CHL.it. X<N«> )  GC  TO  102 

wv=<;aPT<sv:<*7<Na-i)  )-v(N b-d 

ot=2.*ptst/hv 

GO  TO  IjU 

W«=cORT(SV<* (.5*(ET+CT))  >-.5* (BV+CHV) 

jT=?,*nisT/wv 

T  I  Mc  =  T  IfriE  +  OT 

I."  (N0CH«  tO  «2)  GO  TO  3C1 


y  *¥+**■,+  *  V*f¥f¥*** 


1-1 

M=1 

go  ir  n  ji 

i=i*i 

TF(T.PO.N'T)  GO  TC  201 
GO  to  HOI 

GO  TO  (2  j1»  202?2  J3  »20«*,?05»236$  202»2C7»202f  208»<  C^> »  KTYP 
PPFSr.S*  (F(l)tP(Mi) 

GO  TO  2;.T 
0(\ES  =  °(  I> 

GO  tq  2n 

P^^.'JMPd-D+CP) 

GO  TO  ?10 
?rGS-.5» (OF+CP) 

GC  TC-  210 
ppf  <:=;(.  p 
Gr  To  213 
o'-  rc=  pd 

GO  TO  213 

pr£s= . (°p+°{r) ) 

0-0  TO  21  j 
?rt^=.0*frp+p{l) ) 

GO  TO  23  3 

pRlSs.S* (p (7-1) +qp> 

IF(M.£0. (NP-l) i  OFrPRES 

POQT=p&T^. C {PPtS/ICCQ.) *»8XP> +F6K* < 6G V-ACV) 

0t.'*nsO{  Mi 

OSU3F AsAT^n  (G5AG »  ?URN  A  f  N'T 1  >, CRAC) 

P(i')=PirU-PGOT»CT 
IF(3(M),OF.O.)  GC  TC  212 
NOC«=1 

go  To  x:* 

X-TftO-Pt  M)  . 

GUP<-'AiAT<M(G3A0f  BIPKA,NTA?  >1 ,  XPADi 
DCmsAOc.C  COPAO^OSURFA-XRAO^SURFAi  *CRC*GNS) 
CMCS(Nf'-CHAS(N)'OCH 


S2 


GAM/ME/72-2 


r-  UN 


cor  8600  PTN  V3.0-279C  0FT=1 


IF(INtQ.l)  GO  TO  30 
GO  TO  2  A  5 

v  TCHASC=TCM6SC»OM 

XKPaCCMAS/  (CGKA5+CM) 

CGMAS  =  rG‘1ftS«-Dh 

CT-TISO¥(l*-XMk) +CT¥XMR 

VOL  =  CA¥1CHL-X  <L) J *12»- (TCHASC/CRO) 

CHP^-CGMftS/VOL 

CF=12,¥Cr,MAS¥RG*CT/(VCl-CCVOL*CGMAS> 

GC  Tn  ?4 5 

243  TCVASHsTOM AS9-DM 

XMP^OGMA?/  COGMAS^-nM) 

JCMASs“GMAS*Drt 

nT=TTSO* (1 ,*XMR/+PT¥XMR 

VOL'P  A*  (X(L  +  1)-CHU  *12.-nCMA$B/CR0> 

RRG-CCMA5/V0L 

BPrie.^aGHAS^PG^RT/^OL-CCVOLOBGMAS? 

245  GOTO  (250  , ?70, 252,254, 270 ,270 ? 270 , 250,270, 256, 260  ,KTYP 
250  IF(M.£Q.2)  GO  TO  270 

M=? 

?51  OM=OMR 
L  =  I 

GO  TO  230 

252  TF<M.FQ.2)  GG  TO  270 
M  =  2 

O(-:=0MR 
L  =  T 

GO  TO  235 

254  IFIM.EQ.2)  GO  TO  270 

Ms? 

0^-OMR 
GO  Tn  243 

256  IF(M,tn,3>  GO  TC  270 

IF (M.E0.2J  GO  TO  25B 

w-2 

0M=0wC 
GO  TO  240 
25.3  M  =  3 

GO  TO  251 

263  TF(M«E0.3a  GO  TO  270 

TF(MaE0.2)  GO  TO  252 

Mr? 

jP.=0MC 

L”T 

GO  to  235 

262 

OM»OMP 
GO  TO  240 

273  JFtI0.cQ*l>  GO  TO  3t> 

IF(M.CO.  (NP-1H  GC  TO  275 
N-fHi 

GO  TO  im  • 

2?"  IPd.CO.  (U3-1)  ,ANC*KTYP,FG.3.0R«KTYP.£Q*5>  GO  TO  277 
IF ( T ♦ LT» (NS-i>>  GC  TO  280 
GO  TO  30i 


*4 
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COC  6600  FTN  V3.C-279C  0PT  =  i  G2/07/72  1 


‘•'7  nr-rrt'-np 
r."  to  3*u 

28C  T  A V  =  DT 

IF<y  (I)  .LT.CHL.  flNCwCHL.LT. X<I4-1)  >  GC  TO  288 
np=P(I)-OP 
!82  1  =  1*1 

I F ( I . LO . N°)  GO  to  301 

IF (X <I> .LTcCHL. ANO.CHL.LT. XU+IH  GC  TO  292 

134  K=  I  ♦  1 

SV=SQPT(SVK*T<X> ) 

wv=gv+v  a) 

OIST=X(K)~X<I) 

TTC=P1ST/WV 

If-  ( T TC  •  GT •  T A V)  GO  TO  ?01 
T  AV=T  AV-TTC 
P(I)=P(I)+0P 
GO  TO  282 

!33  IF  (KTYP.E0.3.0R.KTYP.E0.5)  GO  TO  290 
OF=OP-OP 
1=1  +  1 
GO  TO  254 

:9c-  op=np-op 

9F=np+QO 
1  =  1  +  1 
GO  TO  ?84 
92  <=I+i 

GV=SCRT<SVK*(,5*(BT+CT)) ) 

MVaSV+t 9*(CHV+By) 

OlSTrXOO-Xd) 

TTC=niST/HV 

IF  (TTC.GT.TAVi  GO  TO  301 

T  AV=T  AV-TT C 

CF=CP+OP 

qr=op+3p 

GO  TO  232 


SHELL  MOTION  SFCTION 


OSHV=?HV 

Pl'cS  =  .5*  (P  (N^-l)  +0P) 

I r  (X ( N1'  - 1 )  .LT.CHL)  P9rS=.5*{9P+OP) 

A0P=PAM1,  +  G AMA  +  CSHV /  (2.  +  SV A)  *  {  (GAMA^l.  >*OSHV/  (2  EVA)  +SCR7  C  (  { GAM 
1 A  + 1 » ) *OSHV/(2.»$VA) )**2.+4.)  )  ) 

$>JA=52.17u+  (PPES-6DP-FP)»nA7SHM 
. =OSHV+SHA+0T 

HX+SHV»SHV/ (SHA*2.) 

/.Oc.GUNL)  GC  TO  1500 
C‘*  )  =SHX 

W(NO,1)=W(N9,1>MSHV-09HV) 

IVJA(NP,1)  =  G 


GAM/ME/72-2 


n UN  CDC  6600  FTN  V2.C-279C  OPT=l  02/07/72 

*  WAVE  PROPOGA1ICN  SFCTION  * 

* 

*?V*«*«***V1M *♦¥♦**¥****#*¥#*¥***  **»<«**»  ¥¥*♦*♦#***#•#<,****<'♦*  *♦♦♦♦»*** 

4vp  r-2 
j=i 

rv=o 

410  IFdt'A  (ItJJ.En.0>  GO  TO  470 

I c ( J. EQ «  2 • OR ♦ J. E9  .4  >  GO  TO  420 
IF (TntO. NO )  GO  TO  430 
1  =  1  +  1 
r,C  TO  410 

4  2ii  IF(T.EO.l)  GO  Tf)  430 
I  =  I  —  i 
GO  TO  410 

43C  GOTO  (440, 450,460*590)  ,  J 
443  J=2 

I=N9-i 
GO  TO  410 
45G  J=3 

1  =  2 

GC  Tn  410 
460  J  =  4 

T  =  N°-1 

r-o  m  410 

47'’  OVs'W(XyJ) 

IF  ( KCHH*  £C| « 1  e  OR » KCHM*  EO •  2)  GO  TO  471 

:*=i 

j^=j 

TAV^T 

471  IF  ( J  •  rO .  2  <  OR « -J«  cO  •  4 )  GO  TC  475 

IP(>:  CI-1)  .LT.CHL./)NO.nHL*LT.X(I))  GC  TO  500 
K  =  T-1 

SV=SOPT(  SVK  +  T  (tf)  } 

WV-5V-V  CO 

Ic(WV»LT.0«)  GO  TC  496 

oi^Tavm-xm 

GO  TO  430 

h75  IF(X(I) .LT.CHL, ANG.CHL.it, X(I+1))  GC  TO  500 
<  =  !  +  ! 

SV  =  SOf'T  <SVK*T(I>  j 
WV=SV+V(T) 

otst=y  co  -xm 

4  3''  TTC  =  nIST/WV 

IF (TTC.OT.TAV)  GC  TO  495 
4fil  I4A(T,J)=l 
M  (I, J}*’, 

TAV=T  AV-TTO 

TF(  J.ECd  .  CR.  J*  FC.3)  1  =  1-1 

CP  =  -,373*R0<T)*S0RT(SVK*T<IH*0V 
°(I) =R(T)+CD 

ovs'im  i  x)  rpvsHf'd)  +nv 

IF  ( J *  E Q «  2 •  OR  •  j .  EO « 4 )  Jssl  +  l 
IHACI,J)=8. 

w(i,j)=wa,j)*iv; 

Trd.EO.l)  GO  TO  405 
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GAM/ME/72-2 


5, 


Tc/(TvC°ANn>  G0  T0  *°0 
(TX.t.O.i)  GO  TO  Cfl2 
r'P  TO  4/1 

<*8?  rv=u 

TWA (T,J) =2 

w  fT»J)=OV 

T-K 

J-JS 

<jO  TO  i,jn 

l  W  A  ( I  f  J )  r  i 

wd,J)Sn. 

Tf-  ( J.  FO.  l)  j-j 
IrU.rQ,3)  J;:4 

W(I,J)=W(ItJ)+CV 
o^ro';?;u  G0  T°  *«* 

86  iy=o 

iwA(r,jjS2 

1*1? 

j=j<; 

Cn  TO  A 1 J 
,0  rwA(i,j)Sl 
wc,  j)  =c, 
t^J.lQ.2)  J  =  i 
IF(J»FQ,4)  js? 

iwfl  a,  j)  =c 

^  T  »  J}  =H  (I  f  JJ  ^.g-y 
r  (IX, £0,1)  G0  70 
GO  T0  47} 

l  rv=n 

XhA(If j,s^ 

T~rs 

j=js 

GO  TO  Ai'5 

Jr  <kTMH.c0.j,  G5  T0 

«  7° 

It«a,IJ)T^T,,,5'0IS7i)  =°  TO  -97 
w d  »J)=OV 

r-is 

j=j? 

r-C  T«  41  n 

ty*t 

GO  TP  4  31 

rv.A  (  T  ,  j|  =2 

f  T*  J)  =0V 
1  =  7° 

J=J'T 

nV  =  OV‘HH{I,j) 

KPHW=? 

- A  1 1  f  J)  zQ 

tav=tac 

rp(i.Fo.NP)  go  Te  493 


COG  6600  FTN  V2.fl.279C  0FT  =  l 


02/07/72 
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GISH 


COC  6600  FTN  V2,C-279C  HPT  =  1  02/07/72 


3 


GO  TO  470 
.c  <CMW=P 

GO  Tn  4 96 

50o  Tf  (KCHW.ro. i)  GO  TO  4' 
rr{KCMW.£0.2)  GC  TO  4' 

rr< j„lQ.2. op.j.eg.^)  c  505 
Ks-l-t 

SV  =  .  5*  (SQR  i  (SV<*CT>* 

WV  =  SV-.GMPV*C,HV) 

Ir(WV .LT.fi.)  GO  TO 
■»lST*y  dj-xi  o 
GO  TO  GO 7 
5GG  <  =  IM 

S  V  =  .  GM'SH^T  (SVK*0T  )  •  T(SVK’ 

WV  =  SV+,G^(8V  +  r.HV) 

IPCWV.LT. ft.)  GO  T f'  •  -(• 

OIST=X (K)-X(II 
Pu7  TTC=nrSf/WV 

IP (TTn.GT. TAV)  GO  >0  496 
t  A\/=  T  AV-TTC 
TfinsTftv 
IKfl (I,J)=1 

L  - 1 

iF(j.ta.i.0P.J.E0.r-i  l=i-i 

9V-nV+0V 

Dr=rH»0*GA  *nV*0\  Me. 
ar.MAb  =  nGMAS  +  riM 
rr.  '!AS=CGmAS"Om 

0‘jP=  co 
OTPrrp 

VOL=  (X  (L*l  )-nwL)¥CA*l?.-*TCMAGR/r 

'5f'  =  i?.*',>G'}ASMG*8T/  (VCL-CCVOL*PG*AS) 

"•■0=bG'iAS/V0L 

(  SP-'OOF)/C’RO*S09T(?VKBBTJ  ) 
VOL* (OHL-X  (L) >*CAM2.~TCMASC/CR0 
CPsi2.*CGMAS*9G*CT/ (VOL-CCVOL*CGHAS> 
C»'-TO  =  CGMAS/VOL 

nvfs-2.68’* (CP-OCM /(CH90*SQ9T<SVK*CT>) 

J  r  ■-  ? 

T  F ( OV G-L T. 9 . )  J9  =  4 
Jf.  =  l 

IFtrvC.LT. 0<)  JG  =  2 
IriJ.t0.i.0R.J.E0.2)  GO  TO  581 
oviUr ( i) =ovsuh( is ♦nvo 
•■>VSUM  (1*1)  =DV5UM  (1  +  3  1  +OVB 
TK-I*i 

CO  TQ  r,62 
581  I°=I 

•TVSUM  ( I )  =  1 VSUM  (  I )  +  PV9 
GViUM (1-1) =OVSUM (I-i>  +  OVC 
I  - 1  -  3.  • 

G .  .  ^TMt4  =  l 

)  =  v;C 

TWA (I,J) =0 
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COC  6600  FTN  V2. £-2790  0F7  =  1  G2/"7/72  1 


M  ( I, J' =W(I , J)  +OVC 
T  AV  =  TAC 
Op  TO  t-7 1 
HO  5^  1  =  1, N13 
no  J  =  l,4 

IF  (TWA  Cl,  J)  .FC.2)  JWA(T,J)=« 


CHANGE.  MOTION  SFCTION 


•**#******  *♦#♦***** 

TF  { NOCH • t  Q  •  2 )  GO  TO  901 
TO  =  ? 

T  =  1 
Ms  i 

no  to  noi 
cin  i=i+i 

IF{T,£0,NT)  GO  TO  601 

nc  to  noi 

A2j  Gn  Tn  in25»630,635,64G,645,S5Qf630»655»630»6CO»£SE),KTYP 
b25  AGV=.5MV(I-1MV<I>> 
no  To  6/o 
6  3"  A  nv  =  \'(I) 

GO  TO  670 

6’5  AGV=«i:#(V(T**i)+CHV) 

nr  T0  670 
'  fiOVa.CO (9V+CHV) 
no  TO  67 n 
64c  Ar-V  =  fHV 

no  T0  670 
ofj  A  o  v  =  p  V 

6C  TO  67 0 

65r»  Ar,V=.5M3V+V<I>  ) 

GO  TO  67" 

£6;  A(,V=.C*<C*V  +  V(T>  > 

00  to  07" 

do5  Anv=,5* (V( r-1) +  nV> 

67:  TF (NCOH.£0.2)  GC  TO  67? 

ACV-CV ( N ) 

IF  (ACV.LF,  0. )  A 6 V  =  C . 
k7?  IF(AGV.Lc. T.)  AGV=0. 

1^(10.60.2)  GO  to  2L0 

no  to  (67^,  691, 70  3,705,720,725, 730,735,740,745,75  0  ,KTYP 
675  PPiF=pci-i)-p(i) 

SGKf's.*!*  (PC  (1-1)  +30(1)  ) 

no  T0  756 

65:  IF(I.cO.i)  GO  TO  690 

TFd.FO.  (N1-l>)  GC  TO  695 
oniF  =  P(i-i)-p<iM) 

6-*n  a<"po=po(I‘ 

nr  TO  7  .>5 

fn  PrTF=P(l)-R<2) 

GO  TP  695 

695  PriF=P(NO-l)-P(I) 

GO  TO  635 
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'HIM  COC  6603  FTN  ¥3.0-2750  OFT  =  l  02/07/72 

pri«*=p  { r-i  i-nP 
AoR0=,5* (RO(I-l)  ♦OHRO) 

CO  t0  79c 

7  *5  I F  ( I  *  i  0.  <MP-  1 M  GC  TO  715 
Pf'IP=.5*  (P(I-l)-p  (I  +  l)) 

710  AGWOr.C* (CRO+CHRO) 

GP  TP  7 55 

715  PuIF=p(T-i)-.c* (RF^CO) 

G('  TP  7  i  1 

72u  PPIF=P(I-1)-«P 
AGkO=PHCO 
GO  TO  75C 

7?5  *>PIF  =  CP-P{I  +  i> 
ftGRO=PRO 
GO  TO  755 

73*  POIF=P (I-l)-CP 
AGRO-PO (I) 

GC  TC  755 
735  PnIF=C°-P(I> 

AGR0=.5MRa(I)  +  <?R0) 

CO  TO  755 

7<*G  PPIF  =  nn-P  ( I  +  l) 

AGRO-PO(I) 

GO  TO  755 
7<o  PO I F  =  C°~ P  ( I ) 

AG’Cs.U*  (CHPO  +  RO  ( I)  ) 

GO  to  755 

75v  pnjFrP (1-1 )-op 

AG,'<0  =  .5T  (PG(T-i)*rPO) 

755  Ai=CMAS(N}/(CRO*CL*12.) 

rV  f  N )  =  C  V'  (M  )  ■*■  (PT  /CMAS(N) )  *  (POTF*61*32. 174  +  »5*AGRC*S1#CC*  (AGV*AGV-2. 

1 *TV (N)*AGV*CV(N) *CV(N> ) ) 

IF(rv(N)  ,(.?«•},)  OV(H) =0* 

Cy<N)sCX(M)+CV(N) *PT 
I  r  ( N .  £.0 «  (N9-1U  GC  TO  801 
\’ ~  N  I 

GC  TP  liil 

* 

♦  GAS  nO  UNO  ARY  »'CTICN  SECTION  * 

*  * 

5.;1  I*? 

305  TF  (X (I) ,LT.CHL.3NC,CHL.LT.XCI+1)}  GO  TO  810 
VfI)=PVSUwiI> 
x  (l>~y (i> *v(I)*ot 
A  =  P  A 

IP  (X  { T  > .LT.CHL) A  =  CA 
r)f'=12.*30tl>  7APV(T)*OT 
GMftc  (T-l)rGMAS(l«i)+OM 

GWAS  (I) =GMAS  il) -CM 
"j7  IF  (T.f  Q.  (NP-1)  )  GO  TO  OOl 
I  -  T  *  1 
GO  T0  *05 
810  QlST=Chl-X<I) 


60 


GAM/ME/72-2 


GI.’M  COC  6600  FTN  V3.G-279C  OFT  =  l  02/07/72 

TTr=pjPT/OHv 

TF (TTC.oT. CT)  GO  TO  320 

TflV-.PT-TTC 

A  =  0GH/>3  +  0G“AS 

X (I) =CHl+nV*OT 

nv?l?M(  I)  ='i V 

OP=12.*CH»C*CA*CHy*0T 

GMAS  (  T)  =  A-OM 

CGMAr--*G“A«  (1-1)  +  9P 

ON=12.*n90*uA+9V*C7 

f(.,MAS=0H 

GMAS  ( I)  =  69  AS  ( I)  -  OF 
r,p  rp 

320  X  d>=X  (I) +CMV+QT 
nvsuwcn =C«V 

T^'  =  12.»  ?0(  T-i)  *CA«CH'/*0T 
G,J  A  S  (T-l)=GMAS(I“i)+OM 
CGMAS=CGMAS-DH 
T  =  I  +  1 

TFd.cn. N*)  GO  TO  901 

X(I)=X(I)+  nV*  GT 

OVS"M(I)=9V 

Ou= 12 • *  °0 ( T) *na#nv»nT 

9  GM  A  c  =  0  GM  A  S  +  O  M 

G“AS  ( I )  =GM  A3  (  H  -Of' 

GO  in  307 

i«**«»**'i*V*«m*«»«»«**444444444444**4444444444444444494«4**444*4444444 
*  4 


*  CHARGE  REDISTRIBUTION  SECTION 


9-'*l  ln-e> 

IF (MOPH.EQ.2)  GO  TO  1000 
L=N°-1 

T  Q  i  C  N=  1 , L 

IC(=MN)  ,GE  .G.)  GO  TO  915 
9 1C  CONTINUE 
N’GCH=? 
go  To  isn 
°15  L  =  N  ° “ 1 

OC  oip  1  =  1, L 
«i6  THAO  ( i) -g  , 

7  C  H  A  r  C  =  9  * 

TCMA9P=G . 
r-i 
‘•=1 

GC  TO  HU 

91«  IFd.FQ.  (NO-1) )  GC  TO  935 
1=1+1 

GO  TO  11 'll 

22.  GO  +  0<o?  <.,922,923,924  /9?5, 926, 922, 923, 922f930, 931), KTYP 
r  '1  7 0‘1A5  ( I-l)  *  TQM A 9  ( I-l)  +CHAS <N> ♦OIST/Cl 

TCHAO  (I)  =TC«ASd)  +CMAS(N>*  (CL-OISTJ/Ci. 

Go  TO  935 

92?  TCMASt I) =TCHAS(I> +CMAP(^> 
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Gun 


COC  6600  F7N  V3.C-279C  CFT  =  1  02/07/72 


GO  TO  o79 

'?  TCMA? (I-i) =TCMAS (1-1) ♦OMflS <N) *01  ST/CL 
TGMAnO=  f  CMASC+OAS  (N>  *  (CL- Cl  ST) /CL 
r,n  to  otq 

<3 2 U  TTMASC-TCM  ASC  +  CMflS  (N)  *niST/CL 

TrMAsn  =  Tc;iASn+CMAS  (N)*  (CL-OIST)/CL 
PC  TO 

9?K  TrMA5C  =  TC‘<  ASC  +  Cf'A*'  (N) 

GC  TO  9 7r» 

9?r'  f  0 M A ° f-=  T »'M  AS n *C M A S  ( N ) 

GO  7"  9 7 9 

92  9  7rMASO  =  TCM/ss9*CKAS(N>  *0 1ST /CL 

TfMAS<I)=TCMAS<I) 4CMAS(N)* (CL-0IST)/CI 
0 O  T0  979 

03  J  TCKASC  =  TC9A$0  +  CMAF(N)*niST/CL 
Tf.MASn  =  TCMAS:T+Cf'AS(N)  *HDIST/CL 
TO,MAC(  I)  =TCMAC;  I)  4CHAS(N)  *  (CL-OIST-WOISTJ/CL 
GO  TP  935 

97 1  TOMAS  (1-1)  sTCMAS  (T-1  )  4CI-AS  (N)*0IST/CL 
TPMASC  =  TCMASC+CMAS(N>  *W0IST/CL 
TCMASn  =  TC9ASn  +  CMAS(N'»  *  (CL-OIST-WOIST)  /CL 

939  I F ( M«  EO •  { N9-1 ) )  GC  TO  1000 
MsM  +  1 

TO  ini 


*  GAS  PROPERTY  REALIGNMENT  SECTION  * 

*  ♦ 


i.j:c  L*Mn-l 

OO  1C7-J  1=1, L 

Ir  (X  (I)  .LT.CHL.ANC.CHL.LT.Xd  +  l) )  GO  TO  102G 
A  =  n  A 

IF(X(  .LT.CHL)  A=GA 

y0L=  (V  (I  +  1)-X(T)  )*A*12.-  (TCMASm/CRO) 

°C ( I ) “GMAS <T) /VOL 

T(I)=o(I) * (V0L-GM0S(T) *COVOL>/(i2c*GMAS(I)*RG) 
GG  TO  i 

1 G 2 C  VOL=  (C.HL-X  (I) )  *CA»12.-(TCNASC/CRO) 
ri-?C  =  CGHA>/VOL 

CT^uF*  (VOL-CGMAS*COVOL)/ ( 12. *CGMAS*RG) 

VOt=  (X  (ItD-CNL)  *9A*12,-(TCMASn/CR0) 

noorl’GvA'v  VOL 

9T  =  nr>* (VOL-RGPAS*COVOL)/ ( 12» *°GM AS*RG > 

1  .•  n  r Cf-!T ? NUr 


PO  T0 


*  CHANGE  LOCATION  SECTION 

* 


¥ 


li'jl  LTYP=C 

Ir(P(N) . G E . C . )  GO  TO  110? 

• ;  GC  w  =  1 

GO  TO  (Iju, 105,615, 620, 91?) ,10 
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CDC  6600  F7N  V2.C-279C  OPT  =  1  32/n-,/7? 


110  5  TF  (I.rn.fJn)  GO  TO  inn 

IF<rx  (N> .GE.X(T) *arO. <CX(K)+CL) .LE.X(I+i>>  GO  TC  1120 
IF'CX(N?  .LT.X(I)  i'U'O.X(I)  «LT#(CX(hH-CL))GO  TO  1115 
1110  MOON* 1 

GC  TO  (100,105, 615, f?0t91*> , TO 
lli^  LvYP=i 
1150  MOrH=f' 

IF  (L  TYP . £0 , 1 )  GO  TO  1136 
IF(T.LE.2.0R*I*G5-  B - 1 > )  GO  TO  1135 
IF(Y(I)  .LT.CML.ft  .-  JIU2)  .LT*CHL>  GC  TO  1145 
IF  {X  ( I -2 ) • GT  #rHL>  GC  TO  1145 
135  I F  ( X  f  I )  .LT.CHL.  ANO.CHL.LT.X(I-H>>  GO  TO  1150 
T  F  ( t  •  E  O  • (N«~l>>  GO  TO  1140 

IF(X<T  +  i). LT.CHL.  AND. CHI, LT.X(H-2>)  GC  TO  1160 
I F ( I « rO. i>  GO  TO  1145 

1140  Ip(Y(I-i).LT.CHL.ANO.CHL.LT.X<I>>  GO  TO  1170 
1145  <TYp=P 

IF(LTYP.EQ.i)KTYP=l 

r. 0  TP  il«r 

1150  TF(LTVp.sn.u)  GC  tq  1155 
KTYP=~ 

IF(CX(M»  .LT.X<I)  .aNO.  (CX  (N)  +CL)  .GT.CHL)  KTYP=11 
GC  T0  1130 
1156  KTYP=4 

IP  (NGCH.cG.2)  GO  ^0  1100 
IF  (TO. £0,4}  GO  TO  113  0 

T  F  ( O  X  { N )  .Gc.X(I)  .flNH.  (CX  (N)  4-CL)  .LE.CHL)  KTYP=5 
IF  (GY  <M)  .G£.CHL.AND.(CX(N)+CU  .LE* X (!♦!>)  KTYP=6 
GO  T0  11  SO 

11*  c  I F  (L  TYP»  EQ , 1 >  GO  TO  1165 
KT yp=7 
GP  Tn  1 1 3 n 
1165  XT Yp=l 

GO  TC  lisr. 

11 7 i  IF(LTYP.FO.l)  GO  TO  1175 
«’TYp  =  9 
GO  TO  1133 
mo  KT YP  =  3 

IFfNCCH. £0.2)  GO  TO  ilflO 

IF(CX(N3 .LT.CHL.ANO. <CX(N) +CL> .GT.X(Ii)  KTY°=1C 
1130  GO  TO  (lj-;  ,62ii,6cJ,62C,213)  ,10 


OUTPUT 


o******«**************«***4.:***««**4«**«*«««««««< 

1351  °fIMT  33  JV  ,TIM£ 

PPIMT  351u,SHX 
D° I  NT  3^20, CHV 
pFJ*(T  35™,SHA 
?°IK'T  354? 

PPJNT  3545* 

L=NG-i 

p-0  1320  1  =  1,  L 

IFti.EG.  (N9-1) ,  ANC.XC  I)  .LT.CHU  GC  TO  1310 
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COC  6600  FTN  V2.C-279C  OFT  =  l  02/07/72  1 


of :P?rc 
V  r  L  =  V  (  T  ) 

ROWrPf>  (  I  ) 

T  =  T(T> 

CSASS=-TCMASCI) 
rr^=v ( I) 
r-n  t n  i32n 

1310  Pf.rc.-np 
np=rp‘s 
v  t  l = "  v 

R0W=rKP 
r  f  m  c  =  r  t 
C,4ASS  =  TCMAS9 
pn«:y (i) 

I3?i)  P^IMT  3590fFOS»cPcStVF-t»TRMPfROWtCMASc 


PETUPN  TO  PRESSURE  PURN  TC  RESTART  CYCLE 


«**<?**¥*» 


GO  TO  lbO 
I5:c  PONTTMUr 


FORMATS 


L  J  J  L 

?  G  *  1 

n  ^ 
«.  y  j  c. 

?j„7 

•j  *> 


FC°M7 T (ftl  , 
FORMAT (h£1 
FO°M  A  T {  12, 
F  C  R  M  A  ^  (2E1 
FORMAT ( 1H1 


3C10 

fopmat ( 1H 

7  1 2  c 

FO  >Mft  T  (  i  H 

V  J  v  w 

POPMAT ( 1H 

FORMAT  (iH 

3  i  5  0 

FORMAT (1H 

-  .53 

FoOH.V'Qrl 

*>-7« 

FOPMAT (1H 

FOPMAT( 1H 

2  22  C 

FORMAT (l“ 

71  G 

rC p*1  A T  (IH 

y  1 1  * 

FM-'MAT  {  !,:■» 

2i2: 

F03UAT  ( ;• .H 

3132 

(  IH 

3  IP  G 

FORMAT ( IH 

3 1 5 1. 

Ff  r!M*T  (IH 

Jlf  0 

F f  -;'**AT  (H 

3 1 7  G 

P02UAT ( IH 

IT*  ) 

T  1 

? C R **  A T  (IH 

3  i  . 

FOPMAT (IH 

32 

c'CH'<At  ( IH 

321  *. 

rO°M A  T ( IH 

*Ain 

S.U) 

I  7) 

•J .  i*  1 

,*GUN  OGSFRIPTTCN*//) 

,*TYPF  CF  GUN* ,T3C  ,  Al 7) 

,*GUN  IFNGTR*,T7C,F13.5,T50,*FT*) 

,*CH AMBER  LENGTF*,T70,F13.5  ,TS0,*FT*) 

,*C.HAMP£F  0 1 AMFT  EP*  ?T70,  FI  3,  F,T«;G,*  IN*) 

,*TPRE  oIAHETER*,T3G,F12.5,T50,*IN*//) 

,*GUN  AST  SHELL  INFORMATION*/) 

,*SWELL  START  pc  ESSURE*  ,  T  2  0  ,F  1 3 . 5  ,  T  SO  ,*  IB  F/SC  U*) 

,  *  GUN  FRICTION  PRESFURE*,T30,F13.5,T50,*LFF/SC  IN*) 

, *SHFLL  MASS*,T3G»F13« 5, T50  ,*LPM*//> 

,*PROFELLANT  INCQPMAT ION*/ ) 

,*TYP£  CF  oROBFLLANT*fT70,  A  1 3  > 

, *  PROF  E  LL  ANT  M ASS* , T30 ,F 13 . S , TSO , *L 0M* ) 
t*  IGNT  T  t  P  MASF*,T3<],FI3.5,TSO,*LOM*) 

,*PPOP£LlANT  DENSITY*, T30,F13. 5, T50  ,*L6H/C»)CIC  IN*) 

, * ISGCHCRIG  FLA  HE  TEHP* , T 3 G ,F 1 3, 5 , T 50 , * CE C  R*> 

,  *FC3CE  CCNSTANT*,T3C-,F13.5  ,T5  0,*FT-LPF/LFH*) 

,*PRESSUPF  PU°N  PATE  CCEF*,T30,F13. E,T5C, *IN/?EC-1000  PS 

,TPCESSL’RF  nURN  RATE  EXFOR ENT*  ,T30 ,  FI 3.  ?) 

, *EROR  I VE  BURN  RATE  COEF* , T30 ,F13. 5 > 
,*CCVOL*JHr*,T3O,F13.5,T5O,*CU0IC  IN/LRM*) 

,*RATIO  OF  SPECIFIC  HEATS* ,  T30,Fi3,  (?) 
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r.'jM  COC  6SQ0  FTN  V3.  C-279C  0PT  =  J  02/G7/72  1 

FCE.MAT(;w  ,*MASS  PFK  G^A  IN  *  ,  T3  0  ,  F 1  3 . 5  »T50  ,  ♦IP?'*  ) 
s  ;i  CCPWH1H  **nRA5  OOEF* , T3G ,E 1 3 . 5// ) 

7?4J  -  CyMAT  ( lu  ,  ’ATMOSPHERIC  ccnm  tions*/ ) 

32?  C  Pf>i}‘iAT(tH  ,  ’PPESFUPF*  ,  T3G  ,Fi  3. 5  ,T5C ,  ’L0F/SQ  IN’  ) 

72*1.  rOR?'Ar(tH  ,*TFMPEFATU°F*»T3G»F13,5,T50,*OFG  R*> 

327C  F C^AT(1H  , ’DENSITY*, T3U,F13.5,T50,’LFM/CUeiC  IN*) 

F C ° M A T  ( 1 H  , ’SONIC  VELOCITY*, T3C , F 1 3. 5 ,T50 , *FT /SEC V /> 
o’9C  FORNATQH  , ’PROBLEM  VARIABLES*/) 

73  j  o  c0?>t£T(in  ,*TIME  INCREMENT*, T3G, FI"*. 5, 750, ’SEC*) 
o?ln  rr?HAT(iH  ,*NUM]CR  OP  GAS  *»PUND A  PIES  * ,T 3 0 , Xi 3) 

•*'Uo  pcp^AT  (1  Hi  ,*PR0CFLI  ANT  GRAIN  9URN  DISTANCE  VS  SURFACE  APEC*//) 

■J3t>C  F  o A  T  ( 1 H  ,*°URN  CIST  (IN)  *,T37, ’SURFACE  AREA  (FQ  IN)*//) 

FORMAT  (IN  ,F13.  5,  TAi«,  F13.  F) 

25.1  eOPIAT ( 1  Hi , *  TIME* ,T?2 , F13 . C , T5C , *SFC*  ) 

35 1  u  FORMA  T ( 1H  , ’SHELL  PCS  I TTCN* * T 32 , F 1 3 . 4 ,T 5 0 , *FT*) 

**52G  FORMAT  (1M  , ’SHELL  VELOCITY*, T32,FiQ»4»T5Q,’FT /SEC’) 

25*0  FORMAT ( 1H  , ’SHFLL  ACOCIFPATION’ , 7 32 ,F10 , 2 , T5C ,* FT /S C  SEC*///) 

■>S'<0  FORMAT  (1H  ,?X,’PO  SI  TION*,T 18,  ’PRESSURE’ ,T33, ’  VELOCITY  *,T4  5,  ’TEMRFP 
iATu=>F*,T6A,*nENSnY*,T77,*PROF  MASS’) 

35U5  FORMA T (1H  , 6X , * ( FT ) * , T 1 5, * (L 9F/SQ  IN) *,T33,* ( FT/SEC )* ,T5G ,* (DEG  R 
1) *,TFof *(L5M/CUF  IN)*,T81,*(LBM) *//) 

355  3  FORMAT  (1H  ,FlG.4,Tl6,F10.2,T3i,F10.2,T46,Fj.Q.2,T61,FlQ.7,T7e,FlC. 

14) 
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COC  6600  FTN  V3.C-279C  0PT  =  1  02/07/72 


ATKM 

R'NCTIC'i  ATKN<X,Y,N,K,XT> 


Amt'N  INTERPOLATING  FUNCTION 


USAGE... 


'rATKNfY ,Y,N,K,XI> 

*-‘>JFR.7..  . 

X  -  TABLE  op  inofpfnofnt  variable  VALUES, 

(HAY  PE  ASCFNQTNG  CR  DESCENDING) * 

Y  -  TABLE  OF  Of FcNOENT  VARIABLE  VALUES. 

N  -  NO.  OF  POINT?  IN  TARLES  X  AND  Y. 

K  -  DEGREE  OF  INTERPOLATION  OESIREC. 

XI-  X-VALUE  FOR  WHICH  INTERPOL  ATI CN  IS  DESIRED  • 

Tuf  TNTEPOLATEC  VALUE  IS  RETURNED  AS  THE  FUNCTION  VALUE. 

Si  CELLS  OF  BLANK  COMMON  ARE  USED. 

DIMENSION  X(N),  Y(N) 

COMMON  II,  Kl,  LI,  LL,  LU 
COMMON  XX ( 13J  *  YY  (13) 

PAT A  KMA  X/  12/ 


IF  (  K  .GT.-  KMAX  ,OP,  K  ,L£.  0  ) 
K1 -K  +1 

T  F  (X(N)-X(IM  10j,10,10 
IP  ( XI-X (1 } )  20,2G,7Q 
LL  =ft 

GO  TP  2L'Q 

TF  (X(N)-XI)  60,40,50 
LL  =  N-Kl 
GO  TO  2C0 
LI  =1 
L 1  ’  =  N 

I1'  (LU-LL-1)  181,180,70 
LI= (LL+LU) /2 
I>r  (XCLIJ-XI)  80 , 80,9n 
L  L  =  L  I 
Gr  TO  ft? 

U'  =  LI 
PC  T0  at 

(XI-X(D)  120,20,23 
I r  CX(N)-XI)  133, 46,40 
11  =  1 
LU-N 

IF  UU-LL-15  183,180,150 
Ll=  (LLHU)*/? 

Ip  (X(LI)-XI)  160,173,170 
IL  =  LT 
GC  TO  163 


GO  TO  300 


ATKNC3C1 
ATKNJ302 
A  T  K  N  3  0  0  5 
ATXNOIC1* 
ATKKjCP5 
AT«n:;C'6 
ATKNCw07 

atkn: jc« 

ATKNC30Q 

ATKKOin 

ATKN-3G11 

ATKN101P 

ATi^NCC13 

ATKN00it» 

ATKNC‘'15 

a^knc^u 

AT^C^i? 
ATKNC01p 
ATKN3?19 
ATK^u^S" 
AT<NCu?l 
AT**- ZZ22 
AT<N0127 

ATKKjOE/. 
A^*\ jZ25 
A  7  K  N  3  n  2  ft 
AT'-'NCTZ? 
ATKN332P 
A  T  <  f.  C  3  2  9 
ATKNi:3r 

atknc: ?1 

Atxa-  j:-32 

a  T  y  \ „  j  3  u 

A?<N£"F6 
A  T  vr  v  ]  r  ■?  7 

AT«n:-78 

atkn;;’^ 

A  t  <  s  ; ;  U 1 
AT < x  ; A. 2 
A  T  X  N  j  :  i.  3 
A  T s  „  -*  i  u 
A  T  X  s 
AT  vs  j<*  Af 
A  t  K  K  j  n  j,  7 
ATKVOf 
ATK\j;t,C5 

AT  XN  :  V5C 
A  T  K  V  G  ”* c  1 
A  T  K  X  c: 5  ? 
ATKN3n53 
AT  2K£C*5<* 
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AT«<n  COC  6600  F7N  V3.G-279C  CFT  =  i  02/07/72  1 


TP 

LL=L  I 

ATXN J055 

0  TO  1  U  ** 

AT<k0966 

13u 

IL  =  LL- (Kl*  1) /? 

ATKS'C997 

IT  f  LL  >  ?0, 21b, 190 

A  TKN1 3  n98 

190 

IF  (LL*K1-N)  200,200,40 

ATKNC359 

?  C  0 

90  21"  I=i,Kl 

ATKNCC6; 

1 1 =  LL  +  I 

ATKKC061 

xy  (i)=x(ii)-yi 

ATKNGf'F2 

21  P 

YY  ( T  )  =  Y  ( 1 1 ) 

A TKN Cl  6? 

90  22u  I  - i  , K 

AT<fNCC64 

90  22 j  J=I,K 

ATN‘.’C0f  6 

22o 

YY<J  +  i>  =  (l#/<XX(J*l>-XX<I})>MYY(I>*XXCJ+i>~YY<  JM) *XX(I) > 

ATKNCrct 

ATKN=YY(Ki ) 

ATKN306? 

Rf  TURN 

ATKNCcea 
A  T  K  N  3  C  c  9 

30D 

oc>INT  10-0,  K 

ATKN3070 

IP  0o 

FORMAT  nH0K=fI12,33H  IS  INCORRECT  FOR  FUNCTION  A  TKN) 

ATKNw 371 

CALL  SYSTEM<2i*0,G) 

ATKN007? 

FWO 

ATKN0073 
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GUN  r-SCRIPTJON 


155MM  MOW 
18. 30000 
‘’.43000 
5.0000  0 
4.18400 


TYCf  or  GUN 
GUN  LENGTH 
fhampep  l“m;th 
CHAMBER  niAMET£» 
BORl  JTAMETcp 


FT 

FT 

IN 

IN 


gun  and  shfll  information 

SHELL  A^T  PRESSURE 
GUN  FRICTION  PRcSSURr 
SHELL  MASS 


4QGG.C0UQ0  LBF/SQ  IN 

350.00000  IPF/SC  IN 

12.770D0  LBH 


PRCFFLLANT  INFORMATION 


TVFE  CF  FRCFt.Lt.ANT 
FRCFELLANT  mas? 

IGM^re  mass 
FRrFELLANT  DrNSI TV 

ISOCHQsjr  'LAMF  TFMP 
FORCE  CONSTANT 
PRESSURE  OURN  RATE  CO.EF 
PRESSURE  SU<?N  RATE  EXPONENT 
ERCblVE  BURN  °ATF  C.QEF 
COVGLU^E 

RATIO  CF  SPECIFIC  HEATS 
MASS  PFR  GRAIN 
DRAG  COEP 


NC  11. jG 
1? » 1600  Q 
.07260 
.05750 
3u00. 00000 
36450G  .0000  7 
.40100 
.67000 
.'jO^ig 
29.62000 
1.4000 
.00214 
.19000 


L  PM 
LBN 

LQM/CU°IC  IN 
9EG  P 

ft-lbf/lrm 
IN/SEC-i QG C  FS 


CUBIC  IM/LRM 
LPn 


ATMCSPHEPir  CONDITIONS 

PRES? UP? 

TEMPERATURE 
0  F  N  S  T  T  Y 

SONIC  VELOCITY 


1 4 . 7  j  GO  Q 
530 .03000 
.  00004 
1128.E5231 


L°F/SO  IN 
OcG  p 

LRM/CURIC  IN 
FT/SEC 


SEC 


FPC9l.EM  VAcIARLES 

TIME  INCRE.'iEmt 

NUMBER  OF  GAS  BOUNDARIES 


00001 
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PROP? LEANT  GRAIN  PURN  DISTANCE  VS  SURFACE  AREA 


BURN  OTST  (IN) 


SURFACE  APEA  (SO  IN) 


•  u  3  ft  •->  c 

•  L.  3  1  5  0 
.  0  2  r-  *5  C 
.  01Q5i** 
.  01370 

.  n12»r 
.omo 
.C1C AC 
.  l'3  90u 

.  ucago 

•  CJ7gn 
.  CGS°C 
. 036GC 
.3J50C 

•  u  C  b  u  n 
.  C  0  ?  1  c 
.00210 
.03110 

0 . 0  3  J  u  0 


1.17700 
1.E67C3 
1 ,3 70  or 

i.A6?oP 

1.5AB3P 

l.r-SCOO 

.87950 

•  7?F7fl 
» 599FG 
« A  BP  90 
.758Fr 
. 3?9SC 
. 27«50 

•  2  ?  A  8  0 
.17810 
.1377^ 

.09130 

•  0 A59C 
0.00000 
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TIMS!  .  C  J 1 7 1  3  SEC 

SHELL  POSITION  2.43''?  FT 

SHELL  VFL'T'ITY  6»6^37  FT/S^C 

SHELL  ACCffLrF1T  I  Cf.  132473, FT/SC  SEC 
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Glossary 

Bore 

Breech 

Chamber 

Chambrage 
Erosive  Burn 

Force  Constant 

Isochoric  Flame 
Temperature 

Muzzle 

Propellant  grain 

Rifling 

Shot  Pressure 

'  • 


The  interior  of  the  gun  barrel.  In  this 
work  the  portion  of  the  barrel  from  the 
area  change  at  the  chamber  to  the  barrel 
exit . 

The  end  of  the  barrel  opposite  from  the 
barrel  exit. 

A  short  length  of  barrel  at  the  breech  end 
witn  a  larger  diameter  than  the  rest  of 
the  barrel. 

A  term  referring  to  the  presence  of  a 
chamber,  as  in  "a  gun  with  chanbrage". 

Propellant  burn  induced  by  the  relative 
velocity  of  gas  past  the  propellant 
surface . 

Term  used  in  the  gun  business  as  a  measure 
of  propellant  energy  potential.  The  force 
constant  is  tne  product  of  the  propellant 
gas  constant  and  the  isochoric  flame 
temperature . 

Temperature  attained  if  a  given  mass  of 
propellant  is  burned  adiabatical ly  in  a 
constant-volume  container. 

The  exit  end  of  the  barrel. 

Small  geometrically-shaped  mass  of  propel¬ 
lant.  A  commonly-used  ^hape  is  c  cylinder 
with  seven  holes  aligned  with  the  axis  of 
rotation  of  the  cylinder. 

A  groove  machined  into  the  bore  to  induce 
a  stabilizing  spin  to  the  projectile. 

An  artificial  pressure  used  in  some  analyses 
(including  this  one).  The  projectile  is 
not  permitted  to  move  until  the  shot 
pressure  is  attained;  this  is  an  approxima¬ 
tion  to  the  force  necessary  to  overcome 
certain  frictional  resistances  to  projectile 
motion. 
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